Quantum chaos and nuclear mass systematics 



O 

o 



in 



(N 

> 

oo 
cn 
o 

00 
O 
m 
o 

> 

I 



X 



Jorge G. Hirsch, Victor Velazquez and Alejandro Frank 
Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, 
Apartado Postal 70-543, 04510 Mexico, D.F., Mexico 
E-mail: frank@nuclecu.unam.mx, hirsch@nuclecu.unam.rax, vicSnuclecu.unam.mx 

The presence of quantum chaos in nuclear mass systematics is analyzed by considering the differ- 
ences between measured and calculated nuclear masses as a time series described by the power law 
1//". While for the liquid droplet model plus shell corrections a quantum chaotic behavior a « 1 
is found, errors in the microscopic mass formula have a « 0.5, closer to white noise. The chaotic 
behavior seems to arise from many body effects not included in the mass formula. 
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The importance of an accurate knowledge of nuclear 
masses to understand diverse processes in nuclear physics 
and astrophysics is well known Though tremendous 
progress has been made in the challenging task of mea- 
suring the mass of exotic nuclei, theory is necessary to 
predict the mass of nuclei very far from stability '2] . Un- 
derstanding of the properties of complex nuclei in terms 
of the actual forces between nucleons is a basic question 
which at present nuclear is unable to answer ^ . For this 
reason nuclear masses are predicted using phenomeno- 
logical models. The finite range droplet model (FRDM) , 
which combines a macroscopic droplet with microscopic 
shell and pairing corrections Q , has become the de facto 
standard for mass formulas. A microscopically inspired 
model was successfully introduced by Duflo and Zuker 
(DZ)U. Along the mean field model it is worth to men- 
tion the powerful Skyrme-Hartree-Fock (HFB) approach 
0. All these mass formulas can calculate and predict 
the masses (and often other properties) of as many as 
8979 nucHdes P|. There is a permanent search for bet- 
ter theoretical models that reduce the difference with the 
experimental masses and produce reliable predictions for 
unstable nuclei. At present, the rms error for 2135 nuclei 
is 674 keV for HFB, 676 keV for FRDM, and 373 keV for 
DZ. The origin of the differences in rms errors between 
the models and the possibility of reducing them are the 
subject of the present investigation. 

Recently, the problem of the mass deviations was ana- 
lyzed from a new angle: in Ref. |^ the errors among ex- 
perimental and calculated masses in were interpreted 
in terms of two types of contributions. The first one was 
associated with a regular part, related to the underlying 
collective dynamics (droplet model) , plus the shell energy 
correction, while the other was assumed to arise from 
some inherent dynamics, possibly higher order interac- 
tions among nucleons (6|, that lead to chaotic behavior. 
According to the latter could be interpreted as remain- 
ing signals of the chaotic dynamics occurring at higher 
energies, whose magnitude suggests that we have already 
achieved (within a factor of 2) the maximum accessible 
precision in the calculation of the masses in mean-field 
theories . It is relevant to ask if this chaotic limit can 
be confirmed by independent techniques and, if so, if this 



lower bound is valid also for mass calculations which ex- 
plicitly include residual correlations, like DZ. 

The presence of chaotic motion in nuclear systems 
has been firmly related with the statistics of high- lying 
energy levels ^ 0. Poisson distributions of normal- 
ized spacings of successive nuclear or atomic excited lev- 
els with the same spin and parity correspond to inte- 
grable classical dynamics, while Wigner's statistics signal 
chaotic motion in the corresponding classical regime [lOj . 
Intermediate situations are more difficult to assess. Very 
recently a proposal has been made to treat the spectral 
fluctuations (5„ as discrete time series Defining 



6n = 



p{E)dE - n, 



(1) 



with p{E) the mean level density which allows the map- 
ping to dimensionless levels with unitary average level 
density, and analyzing the energy fluctuations as a dis- 
crete time series, they found that nuclear power spectra 
behave like 1// noise, postulating that this might be a 
characteristic signature of generic quantum chaotic sys- 
tems. In the present work we implement this idea, using 
the 1/f spectral behavior as a test for the presence of 
chaos in nuclear mass errors. 

In a systematic study of nuclear masses was car- 
ried out using the shell model, in an attempt to clarify 
the nature of the errors. This was achieved by employing 
realistic Hamiltonians with a small random component. 
In jl3l | we have analyzed in detail the error distribution 
for the mass formulas of MoUer et al. ^ and found a con- 
spicuous long range regularity that manifests itself as a 
double peak in the distribution of mass differences [l3| . 
This striking non-Gaussian distribution was found to be 
robust under a variety of criteria. By assuming a sim- 
ple sinusoidal correlation, we could empirically substract 
these correlations and made the average deviation dimin- 
ish by nearly 15% Q. 

In the present letter we carry out a study of the mass 
deviations in the Finite Range Droplet Model (FRDM) 
of MoUer et al. Q , and in the microscopically motivated 
mass formula of DZ 13, Il5l| , analyzing their correlations 
as time series. Two different criteria are employed to 
organize the data, which render similar and consistent 
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power laws. 

To map the mass error data, which depend on the 
charge Z and neutron number A/", in term of variables 
with the maximum possible number of nuclei along each 
chain, the following transformation is employed 

A ^ Int[\/2 (TV sm6l + Z cosd)\ 

Ty_ = Int[\/2 (TV 0056* - Z sine)\. (2) 

Both A and are, by construction, integer numbers. To 
avoid introducing artificial noise, the data are softened by 
the interpolation of mass errors for unphysical values of 
Tj, A, i.e. those with even and A odd, or vice versa. 
This process is necessary to eliminate the large number 
of zeroes which are induced by the transformation, which 
create artificial high frequency noise in the data. 

We found that the best orientation, in order to have as 
many isotopes as possible with the same Tz, is 9 — 56°. 
With this transformation, e.g., there are 174 isotopes 
with Tz — 0. 
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FIG. 1: Mass differences from the FRDM calculations, in 
MeV, as functions of A, for 18 Tz values. 

Figure ^ displays the mass errors AM {A) = Mth — 
Mexp for 18 values of Tz, from = —11 to 6 for the 
FRDM calculations. The regularities seen in Fig. 2 of 
Ref. as regions with the same gray tone are seen 



here in the different plots, as groupings of nuclei with 
similar positive or negative mass differences, for the same 
A region. Besides the two large groups with positive and 
negative mass errors below yl = 50, there are evident 
regions with negative errors close to A = 100, and with 
positive mass differences < ^^O- 
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FIG. 2: Mass differences from the Duflo and Zuker calcula- 
tions, in MeV, as functions of A, for 18 values. 

Figure [21 displays the mass errors for 18 values of Tz 
for the DZ calculations. The deviations are manifestly 
smaller and exhibit considerably less structure. 

The discrete Fourier Fk transforms are calculated as 

where TV is the number of mass differences AM in a given 
series. The parameter 7 makes dimensionless. Given 
that it only affects the global scale of the Fourier am- 
plitudes, we made the simple selection 7 = 1 MeV. The 
Fourier amplitudes are plotted as functions of the loga- 
rithm of the frequency / = k/N for the FRDM data in 
Fig. O and for the Duflo and Zuker data in Fig. ^ using 
a log-log scale. 

These plots have some remarkable features. As ex- 
pected from a l/Zc" power law, low frequencies have the 



3 




FIG. 3: Logarithm of squared amplitudes of the Fourier trans- 
forms of the mass differences, plotted as functions of the log- 
arithm of the frequency, for 18 Tz values, using the MoUer at 
al data. 



FIG. 4: Logarithm of squared amplitudes of the Fourier trans- 
forms of the mass differences, plotted as functions of the log- 
arithm of the frequency, for 18 Tz values, using the Duflo and 
Zuker data. 



larger amplitudes. In most of the plots the largest Fourier 
amplitudes are those with frequencies between 0.3 and 
0.4 (log f = -4 to -3), indicating that oscillations with 
periods A^l ss 20 — 50 are dominant. This is consistent 
with the fitted frequencies found in Ref. A slight 

rise of the amplitudes at the larger frequencies (f « 0.5, 
AA « 2) can be seen in many plots. They represent 
strong fluctuations between some nuclei and their closest 
neighbors. 

The Fourier amplitudes are consistently smaller for the 
DZ data, which have also a Gaussian-like distribution of 
the mass differences |0. This indicates that the FRDM 
mass differences have stronger correlations, which are 
precisely the ones removed in Ref. IT^ . 

The straight lines correspond to the best fitted slopes, 
in the log-log plots, of the power spectra, that is, the 
squared Fourier amplitudes against the frequency. While 
the fluctuations are large, and the number of nuclei in- 
cluded in each chain range from a few dozen to almost two 
hundred, the results are striking and correlated with the 
recently proposed universal features of quantum chaos 




For the 18 chains listed, the slopes are 

"j^izJAf = -1-18 ±0.17, -0.67 ±0.16. (4) 

They fluctuate around -1.2 in the FRDM data and 
around -0.7 for the deviations found by DZ. These slopes 
convey our main result. The former is consistent with a 
frequency dependence of characteristic of quantum 
chaos while the latter suggest a tendency towards a more 
random behavior characteristic of white noise. 

An alternative way to organize the 1654 nuclei with 
measured masses is to order them in a bustrofedon single 
list numbered in increasing order. To avoid jumps, 
we have ordered the nuclei with even A following the in- 
crease in N-Z, and those nuclei with odd A starting from 
the largest value of N-Z, and going on in decreasing or- 
der. Fig. O exhibits the mass differences plotted against 
the order number, from 1 to 1654, taken from Moller et 
al (top) and from DZ (bottom). The presence of strong 
correlations in the Moller at al mass differences is appar- 
ent from the plot. Regions with large positive or negative 
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FIG. 5: Mass differences plotted as an ordered list, taken from 
FRDM (top), and from Dufio and Zuker (bottom). 

errors are clearly seen. In the data of Duflo and Zuker 
the distribution of errors is closer to the horizontal axis, 
and the correlations are less pronounced, although not 
completely absent. 

The ordering provides a single-valued function, whose 
Fourier transform can be calculated. The squared ampli- 
tudes are presented in Fig. El The slopes are 
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FIG. 6: Log-log plot of the squared amplitudes of the Fourier 
transforms of the mass differences, as functions of the order 
parameter(top). Data from FRDM (top) and from Duflo and 
Zucker (bottom). 



"FiiDM = -0-91 ± 0-05, = -0.51 ±0.05, (5) 

for the FRDM and DZ mass differences. 

While this ordering is quite different from the A chains, 
the slopes are very similar. 

To understand the possible origin of these spectral dis- 
tributions, it is worth recalling that, while the FRDM 
calculations involve a liquid droplet model plus mean 
field corrections, including deformed single particle ener- 
gies through the Strutinsky method and pairing the 
DZ calculations depend on the number of valence pro- 
ton and neutron particles and holes, including quadratic 
effects motivated by the microscopic Hamiltonian 
The present results show that the DZ formalism pro- 
duces patterns that are locally smooth approximations 
to the data, and therefore give some information on the 
intrinsic nature of the data fluctuations. 

We arrive at the conclusion that the chaoticity dis- 
cussed in 1^ , according to the criteria put forward in , 
seems indeed to be present in the deviations induced by 
calculations using the MoUer et al. liquid droplet mass 
formula, while it tends to diminish in the microscopi- 
cally motivated calculations of Duflo and Zuker. While 
for the liquid droplet model plus shell corrections a quan- 
tum chaotic behavior a « 1 is found, errors in the mi- 
croscopic mass formula have a « 0.5, closer to white 
noise. Given that both models attempt to describe the 
same set of experimental masses, our analysis suggests 
that quantum fluctuations in the mass differences aris- 
ing from substraction of the regular behavior provided 
by the liquid droplet model plus shell corrections, may 
have their origin in an incomplete consideration of many 
body quantum correlations, which are partially included 
in the calculations of Duflo and Zuker. This interpreta- 
tion would imply that it is in principle possible to reduce 
the limits in accuracy mentioned in [7| for the calcula- 
tion of nuclear masses. It remains to be seen whether a 
robust picture of the coexistence of regular and chaotic 
motion emerge from these studies and whether a quan- 
titative means to evaluate their relative importance can 
be formulated. 
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